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Abstract 

We study genuine CP-odd observables at the LHC to test the CP property of the ZZH interaction for 
a Higgs boson with mass below the threshold to a pair of gauge bosons via the process pp —> ZH 



>^ 



£ + £~ bb. We illustrate the analysis by including a CP-odd ZZH coupling, and show how to extract the CP 
asymmetries in the signal events. After selective kinematical cuts to suppress the SM backgrounds plus an 
optimal Log-likelihood analysis, we find that, with a CP violating coupling b ~ 0.25, a CP asymmetry may 
be established at a 3a (5a) level with an integrated luminosity of about 30 (50) fb _1 at the LHC. 
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I. INTRODUCTION 

The CERN Large Hardron Collider (LHC) will lead us to revolutionary discoveries in particle 
physics. Observing the Higgs boson(s) is arguably the most anticipated discovery at the LHC. 
Once a Higgs boson (H) is observed, a significant part of the LHC program will be focused on 
determining the properties of it. 

One of the most important aspects of the Higgs boson interactions will be its CP property. Due 
to the clear need for new sources of CP violation beyond the Standard Model (SM) to explain 
the baryon asymmetry of the Universe, it is conceivable that an extended Higgs sector may be a 
primary source with or without spontaneous CP violation [1]. There have been significant efforts in 
the literature to explore the possibility to observe the effects of CP-violation in the Higgs sector at 
the LHC [SItJ]- If the Higgs boson is heavy enough to decay to multiple identifiable SM particles, 
then it is quite feasible to construct CP-odd observables. Examples include H — > W + W~, ZZ or 
H — >■ ti |2|. However, if the Higgs boson is light and only decays to light fermion pairs, then there 
will not be enough information for constructing a CP-odd kinematical observable like the triple 
product using the final state momenta of the Higgs decay. The only exceptions are the subdominant 
channel H — > tt, subsequently decaying to at least two charged tracks plus missing neutrinos [3(], 
which suffers from a small branching fraction to the identifiable final state and difficulties for event 
reconstruction. This may lead us to test the Yukawa coupling of ttH. For HVV couplings, one 
may hope to access the clean decay channel H — > ZZ* —> Ai [4] , if the Higgs mass is not too far 
below 2Mz- Attempts have also been made to probe the CP property of the tiH coupling at hadron 
colliders via virtual Higgs effects in ti production [5J , and via the direct associated production of 
ttH at the LHC [6]. However, it is very challenging to establish the signal of the ttH final state at 
the LHC experiments [3]. 

In general, it is known to be quite challenging to construct genuine CP-odd observables at the 
LHC. First, the initial state of the LHC, as a pp collider, is not a CP eigenstate, in contrast to the 
neutrality property of an e + e~ collider or a pp collider. As a result, one may have to seek for neutral 
subprocesses (with qq, gg initial states) or decays of CP eigenstates at the LHC experiments. 
Second, even with the initial states q{pi)q{p2) and g{pi)g(p2), they form a CP eigenstate only with 
special spin and color configurations, and in their center-of-mass (cm.) frame which is in general 
different from the lab frame. These two reference frames are related by a longitudinal boost that is 
unknown and different event by event. Only with large statistics would one expect the difference 
to be averaged out. Furthermore, the symmetric proton beams at the LHC make it impossible to 
identify the direction of a quark versus an anti-quark on an event-by-event basis, which is often 
needed to specify a particle versus an anti-particle. 

In a recent paper [9], some genuine CP-odd observables have been proposed suitable for the 
LHC. One of them is the difference of the scalar transverse momenta [a, |l0|, or equivalently the 
transverse energies, 

0\ = p^, — Pt or E^—E^, (1) 



where p T = \/p^+Py, E T = . /pj, + mi, with superscripts ± specifying the particle charges. 

This observable is CP-odd but T-even 1 and thus generated from CP-violation associated with the 
absorptive part of the amplitude [HI] , thus relying on the existence of an additional CP-conserving 



1 This t is the naive time-reversal transformation which only flips the sign of momenta and spins, without changing 
the initial or final states. 




FIG. 1: Feymann diagram of signal process at the parton level qq — > ZH — > £ + £ bb. 

strong phase sin 5. The other one is a modified triple product 

02 = (Pf x pj) ■ z sgn((p f - pj) ■ z), (2) 

where pc i are the 3-momenta of the particle / and anti-particle f,Pff are their unit vectors, and 
z is the beam direction. This observable is CP-odd and T-odd, and may be generated from the 
dispersive part of the amplitude [llj]. Both variables (H|) and ([2]) are independent of the choice 
of a quark momentum direction and are longitudinally boost invariant, thus adequate for LHC 
experiments. 

In this article, we address the feasibility of discovering CP violation in ZZH interactions using 
these observables. We propose to analyze the process 

pp^ ZH ->■ £ + £- bb with £ = e, \i (3) 

where the Z decays leptonically and the Higgs to bb. The signal diagram is illustrated in Figure [TJ 
The most general Lorentz structure of the Higgs and Z boson coupling, HZ il {k\)Z v {k'i)^ is 

T^(h, k 2 ) = i 2 - [aM 2 z g^ + b (A£fc£ - fci • htgT) + b eT^kiph*] , (4) 

where v = (v2Gir) -1 ' 2 is the vacuum expectation value of the Higgs field, and the Z boson 4- 
momenta are both incoming. The a and b terms are CP-even while b is CP-odd. The simultaneous 
presence of both a (or b) and b in this vertex would generate CP-violation. In the SM at tree level, 
o = 1 and b = b = 0. For a given theory of the SM or beyond, radiative corrections can generate 
nonzero contributions to a, b and b, and thus each of these are energy-dependent form factors 



in general 12}]. In effective field theory language, b and b can be obtained from gauge invariant 



dimension-6 operators such as 



^-H^HV^V^ and -^H^HV^V^ (5) 



respectively [13]. Depending on the nature of the underlying theory, the effective couplings c, c 
could be of order unity for a weakly coupled theory, or of order (47r) 2 in a strongly coupled theory. 
If the theory is valid up to a scale A ~ Anrv, we expect that b and b should be of the order 
l/(47r) 2 to 1. For our phenomenological studies, we will simply take them as constants with some 
optimistic values, and will focus on the complex CP- violating parameter b which accommodates 
both dispersive (real b) and absorptive (imaginary b) CP violation. 

The rest of this paper is organized as follows: In Sec. [Hi we concentrate on the signal and show 
that an asymmetry appears in the observables 0\^ with the above theoretical parameterizations. 
In Sec. IIIIl we consider the backgrounds to this process at the LHC, apply judicial kinematic cuts 
and determine the integrated luminosity necessary to observe these asymmetries. In Sec. IIV} we 
discuss our results further. Finally, in Sec. [V] we conclude. 
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FIG. 2: Total production cross-section for pp — > ZH at the LHC. The left panel contains the cross-section 
versus ruin for b =0 (dotted), 0.15 (dashed) and 0.25 (solid). The right panel contains the cross-section 
versus b for ran = 120 (solid) and 140 (dashed) GeV. 



II. CP ASYMMETRIES IN pp ^ ZH -> l+^bb AT THE LHC 



The CP asymmetries manifest themselves in the distributions of the CP-odd observables Oi } 2 
as introduced in the previous section. We find it convenient to introduce 



AE 7 



El 



E, 



T 



and (pu = sgn((r~ - £~) ■ z) sin -1 ^ x C ■ z) 



(6) 



where <j>u is the azimuthal angle between the two lepton planes £ + — z and t~ — z, multiplied by 
the sign of their logitudinal momenta difference. Our study of the signal is based on a parton-level 
Monte Carlo simulation of pp —> ZH — > £ + £~bb, incorporating the full spin correlations from the 
production to the decay. We use the CTEQ6.1L parton distribution functions 14(. We simulate 
the LHC pp collisions at the cm. energy ^fs = 14 TeV. At lower cm. energies, we would not expect 
any qualitative change for our asymmetry discussions, although the signal cross section and SM 
backgrounds will be different. 

We first present the total cross-sections for the signal versus Higgs mass for 6 = 0, 0.15 and 0.25 
in the left panel of Fig. [2j We see that, for ran = 120 GeV, the SM cross section is of order 700 
fb and b enhances the cross-section by 23% (65%) for b = 0.15 (b = 0.25). In the right panel of 
Fig. [21 we show the signal cross section as a function of b for run = 120 and 140 GeV. There is no 
decay branching fraction included in these figures. We use the illustrative value run = 120 GeV in 
the remainder of this paper. 

To simulate geometrical coverage of the detector and for the purpose of triggering, we impose 
the following basic cuts of transverse momentum, rapidity, and separation on the jets and leptons: 



PTji > 20 GeV, l^yl < 2.5 and ARjjji > 0.4. 



(7) 



With these acceptance cuts, we show the distribution of da/dcpu in Fig. [3] for the case b = 0.25 
(left panel), and the case with b = 0.25i (right panel) by the dashed histograms. The SM result 
is included for comparison by the solid curves in both panels. We see that the CP asymmetry in 
da/d<pn is only induced by CP-violation in the dispersive amplitude as in the left panel with real 
6, but not in the absorptive amplitude as in the right panel with imaginary b. In contrast, the 
distribution of da/dAE? is shown in Fig. U] with the same parameter choice, and the CP asymmetry 
in this variable is only induced by CP-violation in the absorptive amplitude in the right panel with 
imaginary b, but not in the dispersive amplitude in the left panel with real b. 
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FIG. 3: da/d(f>u distributions. The solid curve in both plots is for the SM expectation. The dashed curve 
on the left is for b — 0.25 while the dotted curve on the right is for b — 0.25i. The bin size is 18°. 
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FIG. 4: da/AEx distributions. The solid curve is for the SM expectation in both plots. The dashed curve 
on the left is for b — 0.25, while the dotted curve on the right is for b = 0.25i. The bin size is 10 GeV. 



The asymmetrical distributions in cftn and AEt induced by CP violation leads to nonzero values 
of the cross-section differences 



Act^ = o-0„<o - 0"<fe>o and Aa AEj , = cr A E T <o ~ &AE T >o, 
and in the corresponding asymmetries which are conventionally defined as 



a _ a <t>u<0 -°>ii>0 Q i A _ ff AE T <0 a A£ T >0 

A( t>" = Z TZ and Aae t = Z TZ 



(8) 



(9) 



We show the cross-section difference Ao^ (Acf A e t ) and corresponding asymmetries A^ u (^4a£ t ) 
in Fig. [5] © for a range of b values. The cross-section differences may be sizable and could reach 
about 1 fb for Aer^ and 0.3 fb for Aa A E T with |6| ~ 0.25. The asymmetry in A^ u is typically 
larger than that in A AEt . We see that the cross-section difference follows a linear-dependence on 
|6|, while the asymmetries reach their maxima around Re (b) = 0.25 and Im(6) = 0.25 for A^ and 
Aae t , respectively. This is an indication that the higher order terms in b become significant, and 
thus care needs to be taken for larger values of \b\ when interpreting the results of the asymmetries. 
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FIG. 5: The cross-section difference Aer^,, and the corresponding asymmetry A^, u as functions of Re(6). 
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FIG. 6: The cross-section difference A(jae t and the corresponding asymmetry A/±e t as functions of Im(6). 



III. OBSERVABILITY OF THE CP ASYMMETRIES AT THE LHC 

The signals we are searching for are the ZH events as in Eq. ([3]). We specify these events 
with two jets, with at least one 6-tagged, and two opposite-sign leptons of the same flavor, either 
electron or muon. Contributions to the background mainly come from three sources: 



pp 



QCD bb + 



pp 



tt -»• bb P 



vv. 



pp 



QCD jj + 



(10) 



These backgrounds are calculated using the CalcHEP package [la . Il7 |. The full spin correlation 
have been kept for all pp —> jj£ + £~ where the full 2 — > 4 processes are calculated. In the case of 
the tt background, pp —> tt events are generated and then decayed, ignoring spin correlation. Since 
this process is subdominant, the error resulting from loss of spin correlation should not be large. 
The processes pp — > jj£ + £~ , where j denotes a jet from a light quark (u, d, s or c) or gluon, yield 
large rates due to QCD. Demanding at least one tagged 6-jet substantially reduces this background. 
The tagging and mistagging rates we use are taken from Figure 12.30 of the CMS TDR [18], and 
are listed here in Table [H The mistagging rates are correlated with the tagging efficiency. We 
have taken a rather low value for the 6-tagging efficiency (0.4), in order to substantially lower 
mistagging rates and thus to significantly reduce the backgrounds coming from non-6-quark jets. 



Jet tagging/mistagging rates 



b- quark jet 


0.4 


c-quark jet 


0.03 


Gluon jet 


0.006 


light jet 


0.001 



TABLE I: The probability that a jet will be tagged as a b-quark jet. Taken from the CMS TDR Figure 
12.30. 
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Our signal events should not have much missing energy. We can, therefore, significantly reduce the 
it backgrounds by demanding that the mssing energy in these events satisfies 



i> T < 20 GeV or |^| > 2.5. 



(11) 



We envision the search for this asymmetry occurring after the Higgs has already been discovered 
and its mass is known. We have, therefore, cut the invariant mass of the jets to be near the Higgs 
mass. Further, the b jets from the Higgs decay tend to be energetic, typically around the Jacobean 
peak near Mh/2. For this reason, we apply a higher transverse momentum cut on the jets. Our 
final cuts for this process are 



p Tj > 50 GeV, \Mjj - M H \ < 10 GeV and \M U - M z \ < 15 GeV. 



(12) 



further suppressing the background. With these cuts, the signal over background can be seen in 
Figure [71 We see that the background rate is much greater than the signal for small b. 

The angular distributions of the background processes (see Eq. (110p ) are shown in Fig. [H Many 
different final state partons have been separately presented for clarity. The leading background 
comes from QCD bb £ + £~ , where at least one jet is 6-tagged. The sub-leading background is from 
qg — > qg£ + £~ (labelled as Gji), and the next one is due to ti. The full set of cuts given in Eqs. ([TJ, 
(fTTT) and (fl~2j) are applied along with 6-tagging of at-least one jet. Each background is added to 
the previous one so that the top curve gives the total background. The shapes of the background 
distributions are qualitatively different from that of the signal, while the total rate is much larger 
than the signal for small b (see Figure El) 



p,p->/ + rjj (V7=14TeV) 




FIG. 8: Contributions to the background as in Eq. (jTOJ). Many channels are presented separately as labeled. 
Each background is added to the previous one so that the top curve gives the total background. The full 
set of cuts given in Eqs. (J7J, (fTTj) and (fl2|) is applied as well as 6-tagging on at least one jet. The bin width 
is 18°. 
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FIG. 9: Differential cross section for background plus signal. The solid line in both plots is for b = 
(SM). The dashed line on the left is for b = 0.25 while the dashed line on the right is for b = —0.25. All 
backgrounds are included as in Figure |S] The bin width is 18°. 



Example distributions are shown in the left panel and the right panel of Fig. [9] for b = 0.25 and 
b = —0.25, respectively. The solid curve in each plot is the SM expectation, and the dashed curve 
corresponds to the sum of the signal and backgrounds. It is reassuring to see that the asymmety 
changes sign with a change of the sign of b, which indicates a dominance of the linear contribution 
of b to the asymmetry and thus justifies the leading effective operator approximation. As found 
earlier, the signal asymmetry occurs not far from (fin = 0, while the background tends to be larger 
near (fin = ±ir. This thus motivates us to evaluate the asymmetry A$ u in a restricted region 



126° < (fin < 126° 



(13) 



The resulting asymmetry is shown in Fig. 1101 where all backgrounds have been taken into account. 
We see that the measured asymmetry could reach percentage level. 
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FIG. 10: Asymmetry A$ Ll as a function of b. All backgrounds were included. The inner range of Eq. (fT5|) 
was used to calculate the asymmetry. 



A straightforward Gaussian estimate of the significance is given by 



b ~ 5 A 



\^<t>u\ 



VI 



(14) 



where ot = &4> u <o + a cj>ii>o is the total cross section and L is the integrated luminosity. Inverting 
this gives the integrated luminosity required to discover this asymmetry at an Sa significance 



L = S 2 



ax 



(Acr, 



4>ii' 



\2 ' 



(15) 



The integrated luminosity to measure this asymmetry at la (dotted line), 3a (dashed line) and 5a 
(solid line) level is given in Figure [TTJ For instance, for b = 0.25, one would need an integrated 
luminosity of about 60, 520 or 1440 fb -1 to reach 1, 3 or 5a sensitivity, respectively. 

The sensitivity can be significantly improved if we consider a x 2 distribution of the binned data. 
The "log likelihood" (LL) is defined to be 



LL 



^E 



■/?; 



ndn — + v. 



ri; 



(16) 



where n, is the number of events observed in bin i where f, are expected. The LL determines 
the likelihood that the expected number of events in each bin could fluctuate up to mimic what is 
observed. The larger the LL, the less likely it is a fluctuation. We plot the LL as a function of b 
normalized to units of fb in Fig. [J~2Ta) where we have taken the expectation (i/j) to be that of the 
SM (a = 1, b = b = 0). Once a signal observation is indicated, it would be important to address 
whether one is able to determine the sign of b, which not only checks the consistency of the leading 
operator approximation, but would also give a hint for the underlying physics responsible for the 
CP violating interactions. We plot the LL as a function of b normalized to units of fb in Fig. [J~2Tb) 
but take the observation (m) to be that for a = 1, b = and b = 0.25 while the expectation (z/j) 
is given by the value of b along the horizontal axis. It is interesting to see that the results for 
b = ±0.25 are indeed more similar and it may take more data to distinguish the sign ambiguity for 
a large value of b. 
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FIG. 11: The integrated luminosity required to measure the asymmetry A^ u at lex (dotted line), 3<r (dashed 
line) and 5cr (solid line) using a Gaussian estimate of the significance. 
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FIG. 12: Log likelihood in units of fb as a function of b. The 14 inner bins (from —126° to 126°) were used 
to construct these plots. The left panel shows the LL with respect to the SM for a range of b. The right 
panel shows the LL of b = 0.25 with respect to a range of b. 



In the limit of large rii, the LL approaches a x 2 distribution. To determine the value of the LL 
corresponding with the Gaussian 1, 3 and 5a deviations, we equate the probabilities of fluctuating 
up to that point or higher as in 



Js 



1 



2 / dx- 



-x 2 /2 



dz 



zo(S) 



z n/2-l e -z/2 

2 n/2 r ( n / 2 ) 



(17) 



where n is the number of bins in our LL construction and zq(S) is the LL necessary to give a 
significance of S. If the data is treated with only one bin (n = 1), then the LL approaches the 
significance squared S = V LL. We review the derivation of the LL from a Poisson probability in 
Appendix [Al 

In our analysis, we used 14 bins so that a 1, 3 and 5a deviation corresponds with a LL of 



z (l) = 15.9, 2 (3) = 33.2, zb(5) = 56.0. 



(18) 
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FIG. 13: The integrated luminosity required to measure the absolute value of b on the left panel, and the 
sign of it on the right panel at la (dotted line), 3a (dashed line) and ha (solid line) level. 

Based on these results and matching the LL (fb) in Fig. [T2la) , we can invert to obtain the integrated 
luminosity required to achieve a 1, 3 and 5a deviation of this observable. We present the results in 
Fig- [13| for the integrated luminosity required to measure the absolute value of b on the left panel 
at la (dotted line), 3a (dashed line) and 5a (solid line) level. We find that for b = 0.25, we can 
measure the absolute value of the asymmetry at 3a (5a) with approximately 30 fb -1 (50 fb -1 ). To 
determine the sign of the asymmetry, we must determine the probability that the asymmetry with 
one sign could fluctuate to look like the asymmetry with the opposite sign. Following a similar 
procedure, but using the distribution with b taking the opposite sign for the expectation in each 
bin (the z/j), we present the integrated luminosity to establish the sign of the asymmetry at 1, 3 
and 5<t level in the right panel of Figure [T3l We find that it takes approximately 380 fb -1 (650 
fb -1 ) to determine its sign at 3a (5a) for b = 0.25. 

IV. DISCUSSION 



The ZZH interaction in the context of possible CP violation has been extensively studied for 
e + e~ colliders [20J. Due to the neutrality of the initial state and well-constrained kinematics, it is 
found that an e + e~ linear collider could have significant sensitivity to probe the CP-odd coupling 
of ZZH, especially if beam polarization is achievable. Our work first established the feasibility 
to test the CP violation effects via similar genuine CP-odd variables particularly suitable for the 
LHC for a light Higgs boson well below the ZZ threshold. 

Assuming the observation of a light Higgs boson at the LHC, we have demonstrated the extent 
to which a CP-violating interaction in the ZZH vertex could be explored. It has been shown 
recently that this channel may become one of the Higgs discovery channels due to the improved 
techniques of studying jet substructure |21[] and superstructure [22|, and it is thus conceivable 
that some further enhancement for the signal-to-background ratio may be achieved with more 
sophisticated kinematical considerations. The LHC reach for discovering the CP-violation in Higgs 
couplings reported in this paper should be taken as a conservative estimate. Furthermore, our 
proposal is applicable to any Higgs mass, as long as the production rate is sizable and Higgs decay 
is identifiable. 

Besides the CP studies for a heavy Higgs boson [2|, the coupling VVH for a light Higgs was 
also explored via the weak boson fusion (WBF) production channel [8j. It was found that CP-even 
and CP-odd operators may lead to qualitatively different angular distributions. However, due to 



11 



the lack of particle charge identification, no CP-odd observables can be constructed for the WBF 
processes. 

In this article, we have only studied the effect of b in Eq. (|4]) beyond the tree-level SM. In 
principle, there could also be contributions to a and b which would contribute symmetrically to 
the signal. If this is the case, their fluctuations could mimic small asymmetries requiring a more 
detailed LL analysis. However, the required luminosity for discovery presented in Fig. [IT] would 
still hold. 

V. CONCLUSION 

The need for new CP-violating interactions to explain the observed matter-antimatter asymme- 
try is pressing, and the unexplored Higgs sector may hold the key. If the LHC discovers the Higgs 
boson(s), a detailed study of the properties of the Higgs will then begin. If the Higgs boson is heavy 
enough to decay to W + W~ , ZZ or ti, then it is quite feasible to construct CP-odd observables to 
test properties of the interactions. However, if the Higgs boson is light and mainly decays to light 
fermion pairs, then it would be extremely challenging to test the couplings VVH or Hff. 

We studied two genuine CP-odd observables at the LHC to test the CP property of the ZZH 
interaction for a Higgs boson with mass below the threshold to a pair of gauge bosons via the process 
pp —> ZH — > l+l~ bb. We showed that including a CP-odd ZZH coupling and after selective 
kinematical cuts to suppress the SM backgrounds, we are able to extract the CP-asymmetries in 
the signal events. With a CP violating coupling b ~ 0.25, a CP asymmetry may be established at 
a 3cr (5<r) level with an integrated luminosity of 30 (50) fb _ at the LHC. 
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Appendix A: The Log Likelihood 

The poisson probability of finding n events when v are expected is 

/(n,»0 = -V-. (Al) 

n! 

The likelihood for finding rn events where Vi are expected in N bins is given by 

(A2) 



Locn 



v^e-"* 



To get a normalized likelihood, we divide it by the same form where Ui = rn 



n 



i/.e i 



i>t 



Hi nA l 



The "log likelihood" is then defined as 



rtjln ( — ) + Vi - m 



(A4) 
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In the limit of a large number of events, the LL approaches a \ 2 distribution 

„n/2-l_-z/2 



2 z">- 'e 



X(Z ' n)= 2»/ 2r( n/2) (A5) 

where n is the number of bins N. (If we were fitting to parameters, then n would be N — to.) The 
probability that the expectation value of n would fluctuate up to at least zq or more is given by 

POD 

P(z>z Q ) = / dz X 2 (z,n) (A6) 

In order to determine the zq corresponding to a Gaussian significance of S, we simply set the 
probabilities of fluctuation equal as in 

P(\x-S\ >0) = 2 / dx^=e~ x /2 = P( z >z {S))= dz X 2 {z,n) (A7) 

JS V27T Jz (S) 

and solve for zo(S')- 
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